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Interpolating  Involute  Curves 

Mitsuru  Kuroda  and  Shinji  Mukai 


Abstract.  We  propose  a straightforward  method  for  designing  an  in- 
terpolating involute  curve  whose  radius  of  curvature  is  piecewise  linear  or 
quadratic  with  respect  to  winding  angle.  Designers  can  specify  and  con- 
trol the  curvature  radius  profile  to  a certain  extent.  End  radii  of  a circle 
involute  are  solved  in  terms  of  end  tangent  angles,  and  a G1  involute  curve 
is  derived  by  the  Hermite  interpolation.  For  G2  and  G3  involute  curves, 
relevant  nonlinear  equations  are  solved  by  the  Newton-Raphson  method. 
NC  machines  with  an  involute  generator  can  draw  the  resulting  curves 
with  “reduced  data”. 


§1.  Introduction 

We  present  a new  method  for  designing  two  kinds  of  smooth  interpolating 
curves,  smooth  in  the  sense  of  consisting  of  less  segments  with  continuous 
monotone  curvature  radius  plot.  In  the  method,  we  derive  a G2  (curvature 
continuous)  involute  of  circular  arcs  or  a G 3 involute  of  circle  involute  arcs 
through  describing  its  radius  of  curvature  that  is  piecewise  linear  or  quadratic 
with  respect  to  winding  angle. 

“The  most  important  curve  in  engineering  is  arguably  the  circle  involute 
...  it  has  played  key  historical  roles  in  a variety  of  scientific  and  technological 
applications”  [3].  This  curve  has  excellent  shape  properties  which  make  it 
interesting  for  CAGD  (Computer  Aided  Geometric  Design).  One  can  draw  the 
involute  curve  manually  with  simple  equipment  if  necessary.  NC  (Numerical 
Control)  machines  with  an  involute  generator  are  available  [1,5]. 

In  our  straightforward  design  method,  designers  can  specify  tangents  and 
curvatures  at  junction  points,  and  control  the  curvature  profile  directly  to  a 
certain  extent.  End  radii  of  a circle  involute  are  solved  in  terms  of  the  end 
tangent  angles,  and  so  an  interpolating  G 1 involute  curve  is  derived  span  by 
span  by  the  two-point  Hermite  interpolation.  For  G 2 and  G3  involute  curves, 
continuity  conditions  and  other  requirements  lead  to  a system  of  nonlinear 
equations.  We  solve  this  equation  system  by  the  Newton-Raphson  method, 
using  initial  values  from  the  conventional  G2  cubic  spline  curve.  We  obtain 
examples  of  the  curves  satisfying  additional  requirements,  and  illustrate  the 
properties  of  the  newly  developed  curves. 
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Fig.  1.  Circle  involute  and  its  evolute. 


§2.  Circle  Involute  Arc 

A planar  curve  r(s)  is  expressed  as 


r(s)=Po  + / (l°\ne)ds'  -°0<6,<00'  (!) 

where  s is  arclength  from  the  starting  point  p0  and  9 is  a winding  angle  (the 
angle  between  tangent  vector  and  the  direction  of  the  x axis).  The  following 
relations  hold  among  the  curve  r,  unit  tangent  vector  f,  unit  normal  vector  n 
and  radius  of  curvature  p: 


d2r  dt  d6  1 

ds 2 ds  ds  p 


(2) 


The  radius  of  curvature  of  circle  involute  is  proportional  with  respect  to  6: 


P = 


ds 

d§ 


= Pa+  p{9  - 0o), 


(3) 


P = 


Pb  Pa 

A90 


= const., 


P{0o)=Pa,  p{6l)=Pb, 


where  p is  the  radius  of  circular  arc  that  is  the  evolute  of  r,  and  A is  the 
forward  difference  operator  defined  by  A zt-  = zl+\  — Zi.  We  change  the  variable 
of  the  expression  (1)  from  s to  9 by  the  relation  (3)  and  integrate  it: 


rO  i 

Ap0  = / ptde  = Patio  + p(t  1 - to)  - Pbfii , (4) 

J9o 


r(8o)  = Po.  W)  = to,  n{eQ ) = n0, 
r(^i)=Pu  t(^i)=*i,  n(6»i)  = n1. 

The  vector  equation  (4)  is  understood  easily  as  in  Fig.  1. 
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The  expression  (4)  is  a system  of  equations  with  respect  to  unknowns  pQ  and 
pb-  We  can  solve  this  as  follows: 


_ Aflo  cos(#i  - 4>q)  + sin(fl0  - (ftp)  - sin(6>i  - 4>0) 

^ “ — 2 + 2 cos  A#o  + A0p  sin  A<?o 

_ Aflp  cos(flp  - </> q)  + sin(flp  - (j>0 ) - sin(di  - (j>0) 

^ 6 —2  + 2 cos  A#p  + A^o  sin  A#p  0 

Since  the  involute  arc  is  given  in  terms  of  start  and  end  points  as  well  as  cor- 
responding tangent  vectors  by  the  expression  (5),  we  can  obtain  a G1  involute 
curve  by  the  Hermite  interpolation  which  satisfies  the  equations 


r{0i)=Pi,  i = 0,1, . . . ,n. 


Radii  of  segments  of  the  evolute  of  r are  rewritten  as 


cos(0j  - <t>i)  - C0s(fl;+1  - 4>i)  L 
—2  + 2 cos  A9{  + A (9,  sin  A9, 


i = 0,  1, . ..  ,n  — 1. 


(6) 


§3.  Interpolating  G2  Involute  Curve 

We  can  also  derive  an  interpolating  G2  involute  curve  of  circular  arcs.  Using 
the  equation  (5),  we  can  solve  the  following  nonlinear  equation  system  (7) 
with  respect  to  unknowns  9q,  8i,  ...  ,0n  by  the  Newton- Raphson  method: 


Pi  = P(~9i)  = p(+#i),  * = 1,  2,  . . . , n - 1.  (7) 


However,  the  Jacobian  matrix  necssary  for  the  method  makes  a programming 
code  long  and  convergence  relatively  slow.  Therefore,  from  the  practical  point 
of  view,  we  prefer  to  solve  the  following  equations  (8)  based  on  the  equation  (4) 
directly.  Adding  unknowns  pc,  Pi,  • • • ,pn  to  the  previous  ones  6q,  9\, ...  ,6n, 
we  get 


piTli  + PiAt,  pj_j_i7lj_|_i  — 


i = 0,  1, . . . , n — 1, 


(8) 


Pi  = 


kpi 

Adi’ 


i = 0,  1, . . . ,n  — 1. 


Fig.  2 shows  an  example  of  interpolating  G2  involute  curves  and  its  profile 
of  curvature  and  radius  of  curvature.  Initial  values  were  from  the  conventional 
G2  cubic  splines.  The  Newton-Raphson  method  converged  after  three  itera- 
tions. In  spite  of  the  unpleasant  configuration  of  data  points,  the  curve  derived 
is  quite  smooth.  Its  evolute  curve  (circular  arcs)  is  G1  continuous  except  for 
two  cusp  points  that  correspond  to  extremal  points  of  the  radius  of  curvature. 
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Fig.  2.  Interpolating  G2  involute  curve. 

§4.  Interpolating  G3  Involute  Curve 

Expressing  p as  a quadratic  B-spline  function  of  variable  9,  we  extend  the 
interpolating  G2  involute  curve  in  the  previous  section.  The  radius  p(9 ) is  G1 
continuous.  The  arclength  s is  cubic  with  respect  to  9,  since  p = ds/d9  and 
p{9)  is  quadratic.  In  this  case,  we  get  the  indefinite  integral 


ptd9  — —pn  + pt  + vn, 


dp  d2p 

^ d9 ’ d&  = C°nst ' 

Using  (9)  span  by  span,  we  derive  the  continuity  conditions 
Pini-pi+ini+i-piti+pi+iti+i-i/ini+Vim+i  = Ap;,  i = 0,  l,...,n-l,  (10) 


Pi=p{0i),  Pi  = ~TE  „>  i = 0,  l,...,n, 
av  6=9 i 


d9 2 10=0. 


, i = 0,  1, . . . , n — 1. 
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Based  on  these  conditions,  we  are  going  to  derive  a set  of  equations  with 
unknown  parameters  of  p(9 ) and  solve.  We  use  the  following  notation  in  p—6 
space: 

1)  Knots:  6- 1,  6q, ... , 6n+ 1,  where  the  end  knots  are  of  multiplicity  2. 

9-  i = 0o,  0n+i  — 9n. 

2)  de  Boor  ordinates  [2]:  do,  di,...,  dn+ 1,  where  dt  corresponds  to  the 
Greville  abscissa  (0j-i  + 0,)/2. 

3)  Bezier  ordinates  [2]:  bo,  &i,...,  &2n- 

For  easy  manipulation,  we  break  down  the  non-uniform  B-spline  function 
p{9)  into  the  following  quadratic  Bezier  functions  with  local  parameter  t: 


p(9)  — (1  — t)2b2i  + 2(1  — t)tb2i+\  + t2b2i+2,  (11) 


0<t=^i<l,  * = o,i,...,n-i. 

di+iA0i_!  + diAOi  . 

*>2i  = — — — . * = 0,1, 


A0j_i  + A0j 
foi+l  = di+ 1, 

From  (11),  we  obtain 


i — 0,  1, ...  ,n  - 1. 


_ , _ 2A&2,  •_  n 1 

pi  — 02ij  Pi  — > 2 — U,  1,  . . . , 71, 


_ 2A262i 

Vi  ” (A0;)2  ’ 


i = 0,  1, . . . ,n  — 1. 


We  solve  the  equations  (10)  with  unknowns  9q,  9\, . . . , 9n,do,di, . . . ,d„+i  by 
the  Newton-Raphson  method,  using  initial  values  from  the  conventional  C2 
cubic  splines.  The  number  of  equations  is  2n,  while  the  number  of  unknowns 
is  2n  + 3.  Accordingly  we  can  give  3 more  additional  requirements.  Since  the 
radius  of  curvature  p(9)  determines  a unique  curve  shape,  we  can  specify  and 
control  an  interpolating  curve  by  the  control  polygon  (Greville  abscissae  and 
de  Boor  ordinates)  of  p{9).  Therefore  the  curve  includes  circular  arc,  circle 
involute  and  involute  of  the  circle  involute  because  p(9)  is  a quadratic  B-spline 
function. 

Fig.  3 shows  an  example  of  an  interpolating  G3  involute  curve  with  the 
same  data  points  and  the  same  end  tangents  as  in  Fig.  2.  The  computation 
converged  after  four  iterations.  The  evolute  of  this  curve  has  cusp  points 
within  segments,  since  the  radius  of  curvature  is  quadratic,  while  the  evolute 
in  Fig.  2 has  cusp  points  only  at  junction  points.  The  evolute  of  the  evo- 
lute (circular  arcs)  has  three  cusp  points.  The  curvature  profile  shows  the 
smoothness  of  this  G3  involute  curve. 
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§5.  Some  More  Numerical  Results 

To  illustrate  the  properties  of  the  newly  developed  curve,  we  show  some  more 
examples  of  the  curves.  The  same  data  points  are  used  in  Figs.  1 to  5. 

The  G2  involute  curve  is  practically  more  important  than  the  G1  and 
G3  ones.  Accordingly,  in  Fig.  4 the  G2  involute  in  Fig.  2 is  compared  with 
other  curves:  (a)  the  G3  involute  curve  with  the  same  end  tangents  in  Fig.  3, 
(b)  the  conventional  G2  cubic  spline  curve  which  is  used  as  an  initial  curve 
by  the  Newton-Raphson  method,  and  (c)  a G1  biarc  curve  derived  by  mini- 
mum difference  between  curvatures  of  two  arcs  [4].  The  labels  “G2”,”G3"  or 
“G2”  in  Fig.  4 point  out  which  side  the  corresponding  curve  passes  through. 
Small  circles  in  Fig.  4(c)  are  centers  of  circular  arcs.  It  is  understood  from 
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§6.  Concluding  Remarks 

We  proposed  the  design  method  of  up  to  G3  interpolating  curve  as  an  involute 
of  circular  arcs  or  an  involute  of  the  circle  involute  arcs.  This  straightforward 
approach  provides  a tool  for  the  construction  of  planar  curves  consisting  of 
segments  with  monotone  curvature  radius  plots  of  constant  sign.  Available 
NC  machines  with  an  involute  generator  are  able  to  draw  the  objective  curves 
with  reduced  data. 
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